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Приблизний анал!з неллеийно! 
конвективно! математично! модели 


Изучаются теплофизические процессы, сопровождающиеся фазовыми переходами вещества, описываемые 
математической моделью, в которой температура каждой из фаз удовлетворяет уравнению переноса тепла 
со своими теплофизическими коэффициентами, на границе раздела фаз, обе температуры постоянны и 
равны температуре фазового перехода, а на заданных частях границы поддерживается определенный 
режим. Поверхность раздела фаз («свободная граница») является неизвестной и для ее определения 
дополнительно задается условие Стефана. Это условие превращает математическую модель в нелинейную 
проблему большой трудности. Для описания поля скоростей в жидкой фазе используется система 
уравнений Навье-Стокса. Для решения задачи предложен метод малого параметра. 

Ключевые слова: дифференциальное уравнение, свободная граница, численные методы, 
функционал, оптимизация. 
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Вивчаються теплоф1зичн! процеси, як! супроводжуються фазовими переходами речовин, як! описуються 
математичною моделлю, в який температура кожно! фази задовольняе р1внянню переносу тепла 31 сво?ми 
теплофтзичними коефипентами, на границ! розподллу фаз, температури постйн! 1 дорвнюють температур! 
фазового перехода, а на заданих частинах границь шдтримуеться заданий режим. Поверхня роздлу фаз 
(«в1льна границя») е невдома 1 для Й визначення задаеться умова Стефана. Ця умова перетворюе 
математичну модель в нелйну проблему велико! тяжкости. Для описання поля швидкостей в р1диннй фаз! 
використовуеться система р1внянь Нав’с-Стокса. Для розв’язання задач! запропонован метод малого 
параметру. 

Ключов! слова: диференщальне равняння, в1льна межа, чисельн! методи, функщюонал, оптим1заця 
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1. Принято считать, что чем ближе модель к действительности, тем точнее прогнозы 
и тем эффективнее, следовательно, управление. Однако это не так. Процессы ЭШПИ 
(электрошлаковый переплав) настолько сложны, что, попытавшись построить мате- 
матическую модель, весьма близкую к реальному процессу со всеми его деталями и 
особенностями, можно прийти к очень сложным уравнениям, решения которых крайне 
затруднительны и приводят к существенным ошибкам. Исходя из этих соображений, 
необходимо стремиться к построению сравнительно простой математической модели 
процесса ЭШП, отражающей его самые существенные стороны. 

В настоящее время известны основные качественные зависимости протекания 
процесса ЭШП. К сожалению, полных математических моделей его не существует. 

Перейдем к построению математической модели процесса кристаллизации. 


Будем обозначать через ©“, область, занятую жидкой (твердой) фазой в момент 
времени г. При этом ©, - заданная область в К”, граница которой состоит из двух замк- 
нутых связанных поверхностей Г’ и Г, ‚не имеющих самопересечений, где Г,- гладкая 
замкнутая поверхность, лежащая внутри ©), ‚ такая, что Г, лежит внутри ограниченной 
области, границей которой является Г,. Требуется определить области ©", и ©’, (т.е. 
границы Г’ и Г’), вектор скорости Ус, 0 = (У: (х, Г), У, (х, Г), Уз (х, Г)), давление 


РР), концентрации примеси с(х,!), температуру жидкой и’(х,г) и твердой и (х,г) 
фаз по следующим условиям: 


ВИ м р + (УУ)и* (х, 1) а Уи (х, 1) =0, (х,)ЕБ,, 
А Ума) =0, (х,) ЕД, , 
о +(УУ Ух.) +Ур(х, п) = УУ?У(х, п) + Г(м",с), (хде В", = (1) 


УУ(х, 1) =0, (хр еЕД,*, УУ(х,0) = С); 


7’; 


ТУ, ри =-а(х,п, (хр ЕГ-:У =-а-Р 


+ 


У =0, (х,0ЕГ,, и*(х, 0) = В*(х,), (хрЕГЛ Г, ; 


и* (х,0) = А*(х):и* =и =Т'-5с, Ё О НН ХР", Х,0ЕГ,, 
- бп дп 
О) + (УУ)с(х, 1) -УУс(х, В) =0, (хр ЕР); 
с(х,0) = 8,(х), с(х,1) = 8(х,0, (хр ЕГ), 2 = Вс, (хр еГ,. 
т 


Здесь х=(х,х,,х,), О» ={(х,1): хЕО*,,КО,Т)}, О*, — области соответственно 
жидкой и твердой фаз, 00" =Г, Г’, 907 =Гу чГ,, п — нормаль к Г, ‚ направленная 
в сторону О. Далее, У = (9/дх,0/0дх,,0/0х,), Т(У, р) — тензор напряжений с элемен- 
тами Т, =-6,Р-+у(О, / дх, +0И, /дх,) ‚ У, и У. — нормальная и тангенциальная состав- 


ляющие у я У’ - скорость движения онта кристаллизации в направлении нормали п, 
п 
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Ту, Х.р’,р ,а,В,у,К.,К_ — известные положительные постоянные. Отметим также, 

что если Ф(х,) =и“(х-+Ес(х,Р)-Т’ =0 — уравнение поверхности Г, ‚ тогда \, =-Ф, |УФ|. 
В дальнейшем удобно условие Стефана представить в следующем виде: 

Кии Г.) = К° | Уи Р-К, * | Уи" + (К? + КК, (Уи „Ус -= (+ КК, (Уи „Ус + 

+Хр* (Ки, +К и) + хр"= (К, +К)с, =0, (х,ЕГ,. 


2-9 


В работе предполагается, что А(х)е Н“"® (©), Со) ЕН” (0), В*(хбЕН — 


2+“ 


2 (Г*х[о,Т]), «0 еН“ (©, ). 


(ГА Гу х[0,Т]), Р(и,с) Е СЕ”), (хрен 
1+4 

При этом 8(х,1) и в, (х,0 должны быть функциями класса Н "2 (Взх[0,Т]). 
Считается также, что выполнены условия согласования до первого порядка вклю- 
чительно, которые формулируются аналогично [1, с. 268, с. 363]. 

Цель работы -— исследование задачи Стефана с учетом конвекции и примесей в 
жидкой фазе. 

Постановка задачи — построение приближенного решения задачи, используя 
метод малого параметра. 


2. Будем искать свободные границы Г, и И. в следующем виде: Г, ={х= 
‘= х(@)+и()р(®,0)}, Г’ ={х=х@0)+1@,07@)}, где © = (®,0,),0 = (@,,0,), х(®) ЕГ., 
х0)=Г,’, р(®,) и 1(0,1) -— некоторые функции соответственно классов 
2+9, 2-9, 
и (Г х[о,Т]) и ах (Го х[0,Т]), р(@,0)=0 ип@,0) =0. Введем также 
обозначения О„^ =О, х[ о,Т], Гу =Гь х[ О,Т], Г =Гу ж[,Т], Га =Гох[0,Т]. 
Далее, для достаточно малых чисел = будем искать решение задачи (1) в виде 
следующих разложений: 


и* (хе) = ис (х)+ УЕ и,* (х,Р) ров) = Ро(х)+ У =“ р, (х, Г), 


КЕ К=1 


У (я, =У (+ У =“, (ХВ, сх, = 6,0) + у ве, (х,Р) (2) 


К=1 КЕ 
1=1,2,3; р(®,Е=)= ое р, (®,г), п(@,Е:=)= = ц, @,2. 
К=1 
В работах [2-7] изучены нулевые и первые приближения задачи (1) для малых 
чисел &. При этом установлено, что и,*(х) = А*(х), Уо(®=С(®, в: 2 (х), 
рае + а На, —- + 
р. (®,еЕН 2 (Гу), ПОЮ еГу, и(БрядЕН' ово р ЕН" 7 2 (©>), 
причем р и 1) находим как неподвижную точку сжимающегося оператора М, 
ди _ к 9 
_ дп “, дп 
О ие место следующие формулы: 


+ Х(х,г))аг› ХФ) ЕГу. 


и - 


ди ди ди 
И г, = их += (Л +Ви +—Ю +=” (а, Л, + Вьм, +2) +...+5 “(о +В, +) +0(Е), 
Ох Ох Ох 
(хр) ЕГ а ро А РОО. 
а М” 
(хх еЕГх, 
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Кии ГГ, в) | = К? [Ушо- [ -К,* | Ум" 2] +8[2к_?(Ущ ‚Ум, )- 2%, (Ушу ‚Уи )+Ф, + 
хр: беж. Е, ЕЕ [Уи Уи ,)- 28 (Уи, ‚Уи’,)+Ф,+ур' (Ки, +Ки, )] + 
+0(5“)=0, (хЮеГу, где а, (х,г), В, (х,г), Р.(х,1), Е,(х,Ю) и Ф,(х,Г) — известные 
функции. Из последней формулы следует, что К? |Уи, Р-К,* |Уи |= 0, хЕГ,, 
ди ый ди,” 
- бп * дп 
4. Введем обозначения: М(Е,С,)=<&УХ, КЕУБ, ‚+.-НЕМС,, МЕ, )=ЕМЕ НЕМА, +6 МЕ, 
МУ, рт=Т(,, ртр... рут. 


[4 +Ф, = яр” 9РЬ, (О еГь. 


Затем рассмотрим А-е приближение (У,, и“, р,,р,.П,,с,) задачи (1) для малых 
чисел 5. Имеем: 


= +М,(У,У,)+Ур, = ата? Аш", (хВЕО,* 
1 { 


УУ, = 0, ИО, ; М(У.р,)п = 0, (хЕГ ь (3) 


ОГ? 


+ 


У, (х,0) =0, У, =а-Р — 


Уи + + 2, (х, 0], У, =0, (х, Е Гог 
р’ |Ущ | 


Ам, иг) = а’ Уи, *, (х, 1) Е Ог', 
Ра узи, = 0, (х,Е) Е О | 
1 


оао нию ое ор сие 


| Ум (х(®)) | р, (®, В +и, (х(), Е) + 1, (6), 0) = 0, (х,)ЕГу 


Сем, 9, ср = 0, (х, Е) Е О’, с(х,0) = 0, ск (х, 1) =, 


(ХЕГр; -а о 5) в +Е, ,(х,ЕГур, (5) 
дп [ Уи, | 

дс р 

5 Пк + с, @),1) + в, 098 ет 


здесь Е, (х,Г), р (х,1) и Е (х, Г) — известные функции. 


Зададим теперь У=У, (х,0). Затем решим задачу (4), (5) и найдем и’, с, р, 1.. 
После чего решим задачу (3), являющуюся начально-краевой задачей для системы Навье- 
Стокса. Затем, используя новое значение У, (х,1), снова решаем задачу (4) и (5) и т.д. 


Следовательно, получим процесс последовательных приближений У,,и,”, с,, р,, И». 


Доказательство сходимости этого процесса аналогично приведенному в [8] причем при 


2+, 2 


2+0, 


Е > т 2-9, ‚ ы + тю СИР 
заданном р, (©)еН  ' (Г) найдем функции и,*(х,1,р,)ЕН `? (0), <с(5Ьр)ЕН `? (0) 
как единственное решение задачи (4)- (5), ар, (©,Г) находим как неподвижную точку 
сжимающегося оператора М, : 
ди, _ р ди,” 
дп * дп 


1 1 
М.Р, = +Ф, (х, 0) 4, (х(6), 1) Е Ге. 
0 
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Справедливы следующие утверждения. 
д 
Лемма 1. Пусть выполнены условия | УА* (х(6)) |= — ‚ ХЕ 
и 


некоторая постоянная [1, с. 364] и пусть У’А*(х)=0, хе о’ ‚ А"(х)| = В*(х,0), 


хеГо* Го. 
АСЕ =0, С) =0, хебо, К | УАТС = Е, |УА* О, |УАОЕ>=>0 


(здесь 5 — некоторая положительная постоянная). Тогда оператор М > действующий 
2+4 


изн“ > (Г) в Н Е (Гг) ‚ имеет там неподвижную точку. 

Лемма 2. В качестве К-го приближения задачи (1) можно взять решение 
и“ (х,1), с,(х,1), У (х, 0), р, (©,1), |, (@ г) задачи (4) — (5). 

Теорема. Пусть выполнены условия леммы 1. Тогда приближения У, (х,0), и, (х,0), 


с,(х,1), р, (6,1), 1, (© ,г) сходятся к функциям У(х,г), и*(х,г), с(х,г), р(®,1), п, (@,1) — 


2+а,—— 


класса Н * ‚ являющимся решением задачи (1). 


де, (х) 26 
[2 


Лемма 3. Пусть на Гу’. Тогда при малых числах = и достаточно 


малых значениях {справедливы формулы: 


5 с. (х(@9),1)- 8,(х(0),1) 
0 
дп 
т 5%: (0), Г) + № (@),1) 
Г: х=х(о) пре ОИ ОТ 


Замечание. Доказанная теорема фактически устанавливает существование реше- 
2-9 


г: [4 
Г” :х=х(0)-п»`= 


1=1 


+0(5^), ера 


+0(5*), хеГ.. 


ния задачи (1) в классе функций Н "72. 
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ОУ ло 


А.5. МтепКо 
Арртомтапоп Апа[у515 о7 Мотптеаг Майетанса! Моае 
ший СопуесНоп 
ТЬ1$ рарег ежеп4 ю Фе ите-4дереп4епЕ сазе зоте гезий оМаштеа Бу Ше аиог$ Гог 


(Бе уеаду-попз(а{е ЗерВап рго Мет \ сопуесйоп. КеГегиие тю Ше 1се-\уаег зузет, 
\аег 1$ аззитеа {о Бе шсотргезЫе апа {ю оБеу Ше ЗюсК$ едиайоп: 


(хо: 12 _ 
Е. (УУ)У (х, 8) + Ур(х, в) = ИСИ) Ра 
01 Бе 
(Ув уеосипу, У 13 Кшетабс у1зсозИу, р 15 ргеззиге, {15 Биоуапсу огсе, и" 15 \аег 


{етрегавхге), \Б|е {етрегаге и’ занзНез Ше Беаё сопаисНоп-сопуесНоп едиайоп Уи 
{етрегашге-4ерепаеп( Фегитрруз1са1 ргорегиез: 


+ 
— + 
и (УУ)и* (х,1) -а2У?и*(х0=0, (х0 ЕТ 
Тре (етрегайте йе! и шт Ше зоПа рВазе 15 соуегпед Бу Ч из1оп ощу: 
а. а2У^и` (х,г) =0 (х,0Е ОТ. 
01 
Ас Ще 1се-ужег пщегЁасе, и* =и =0 апа МеРап сопд1юоп Во148: 
и = 0, У К. — К. и со$( и, х,) + К с0$( п,Е) = 0, хеГ,, 
= Ох, Ох, 


+ 
р = (х.0 а = а. Е (0,1) О ФЕ 

ТБе зсВете 1$ сотр!ееа Бу па апа Боипдагу соп 11 оп$. 

Тве аиог ргезеп5 ап ех1${епсе Феогет Гог Ше сазе оЁ гее зрасе 4пптепз1опз. Тве 
таш АЁЯсоу соп$1565 ш Ше Гасё Фа ю пиегрге Фе ЗюКез едиайоп ш а \ееК зепзе, 
зоте шГогтайоп 15$ пееде4 оп Фе геглоп у\/Беге Фе {етрегаеге 15 розйуе, \мсь 1$ ш 
(иго шЯчцепсед Бу Фе уе]осиу Не!4 изе!. ТВе ргес1зе Гогти]айоп оЁ Фе рго ет гедиге$ 
а (есбса| сВо1се оп Рапсйоп $расез. Ех1$епсе оЁ а зоаНоп 15$ ргоуе4 Бу шиодисше а 
{етрегашге 4ерепдепе репаКу {егт ш фе Нша По\ едиайоп ш от4ег {ю Аейпе Бой Фе 


арргохипайпе ‘етрегавите И * ап@ Фе арргохипайоп уеюосйу У ш#е мВоЕ дота. 
Сотрастез$ агхитеп! аге изе4 {0 оеё а сопуегоепЕ зибзедиепсе, \по]е Ши 1$ зВо\уп ю 
зоуе Фе опетша! ргоет. ТВе диезНоп оЁ ип1лацепез$ 1$ [ей ореп. 


Статья поступила в редакцию 02.12.2011. 
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